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ABSTRACT 


To  study  the  rotational  motion  of  projectiles  in  free  flight,  yaw 
sonde  devices  are  mounted  in  them  to  record  the  motion.  Since  these 
devices  use  the  sun’s  rays  to  obtain  a  record,  it  is  necessary  to  fire 
test  rounds  only  at  times  when  the  sun  is  in  a  favorable  position. 

This  report  describes  a  method  for  determining  these  times  for  any 
given  firing  site.  Year  charts  for  three  sites  (NWC,  China  Lake, 
California,  Yuma  Proving  Ground,  Yuma,  Arizona;  and  NASA,  Wallops 
Island,  Virginia)  are  shown.  From  such  charts,  the  most  favorable 
firing  times  may  be  chosen  for  each  day  in  any  given  year.  Range 
azimuth  is  the  most  important  factor  in  determining  whether  the  yaw 
sonde  will  see  the  sun.  The  nearly  north-south  range  at  China  Lake, 
California,  is  most  favorable,  allowing  firings  before  0830  hours  and 
after  1430  hours  all  year  long.  The  east -west  range  at  Yuma  Prr  ing 
Ground,  Arizona,  is  least  favorable  allowing  firings  only  in  three  win¬ 
ter  months  between  0900  and  1400  hours.  Wallops  Island,  Virginia  al¬ 
lows  firings  after  1300  hours  from  September  to  April.  Out  of  the  0* 
30*  ,  45*,  and  60*quadrant  elevations  considered,  the  60*  quadrant  ele¬ 
vation  offers  the  most  restricted  opportunities  for  firing  at  every  site. 
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NOTATION 


A 

Sun's  azimuth 

E 

Sun's  elevation 

D 

Sun's  declination 

L 

Latitude  of  firing  site 

n 

Range  azimuth 

Q 

Quadrant  elevation 

s 

Unit  sun  vector 

V 

Unit  velocity  vector 

f 

Angle  between  S  and  V 

LHA, 

H  Hour  angle  of  sun 

CVT 

Civil  standard  time 

LMT 

Local  mean  time 

DLG 

Longitude  correction  on  LR< 

MP 

Meridian  passage  at  LMT 

V 

Initial  (muzzle)  velocity 

X 

Horizontal  coordinate 

y 

Vertical  coordinate 

g 

Acceleration  due  to  gravity 

r 

Horizontal  range 

h 

Height 

y' 

dy/dx  =  slope. 

Sub  sc 

ripts 

1,2,  3  On  any  vector  denote  its  x,  y,  z  components 
X,  Y,  Z  Coordinate  axes  in  Figure  I 


? 


INTRODUCTION 


One  of  the  devices  used  to  study  the  angular  motion  of  projectiles 
in  free  flight  is  the  yaw  sonde  (Ref  1).  Since  these  devices  use  the 
sun's  rays  to  obtain  a  record,  test  projectiles  should  be  fired  only 
when  the  sun  shines  from  a  particular  direction  on  the  side  of  the  pro¬ 
jectile  while  in  flight.  To  schedule  yaw  sonde  firings,  it  is  therefore 
necessary  to  know  in  advance  when  the  sun  will  be  in  a  favorable  po¬ 
sition.  The  work  described  in  this  report  was  done  to  fulfill  this  need 
at  three  firing  sites.  A  general  method  was  devised  for  predicting 
the  sun's  visibility  at  any  given  range  (azimuth  and  latitude)  on  any 
day  of  any  year. 

The  computational  techniques  employed  here  are  more  compre¬ 
hensive  than  those  used  previously.  La  Combe's  technique  (Ref  2) 
was  restricted  to  a  single  day  in  a  certain  year,  based  on  Nautical 
Almanac  data  for  the  date  in  question.  The  basis  for  our  new  tech¬ 
nique  was  furnished  by  Doan  and  Sandford  (Ref  3)  who  showed  that  for 
many  practical  applications,  where  an  error  of  1*  in  the  sun's  position 
is  acceptable,  the  year-to-year  changes  in  the  sun's  position  on  a  given 
day  are  negligible.  Doan  and  Sandford  provide  two  sets  of  graphs  the 
combination  of  which  forms  the  key  ter  the  solution  of  our  probletW/  m. 

Two  independent  lines  of  computation  are  combined  to  obtain  the 
desired  results.  The  first  line  uses  the  yaw  sonde  trajectory  charac¬ 
teristics  to  determine  where  the  sun  must  be  for  yaw  sonde  record¬ 
ings.  The  second  line  uses  the  chart  devised  on  the  basis  of  Reference 
3  to  tell  when  the  sun  will  be  in  these  desired  positions. 


THEORY  AND  PROCEDURE 


Limit  Calculations 


The  purpose  in  this  first  line  of  computation  is  to  locate  the  sun 
when  it  is  on  the  edge  of  the  field  of  view  of  the  yaw  sonde  at  certain 
stages  in  the  flight  of  a  projectile.  The  position  of  the  sun  is  defined 
by  a  unit  vector  S  which  is  defined  by  two  angles,  A  and  E,  represent¬ 
ing  azimuth  and  elevation,  respectively,  as  Figure  1  shows.  In  the 

same  figure,  it  can  be  seen  that  the  components  of  this  unit  sun  vector 
arc 
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cos  E  cos  A 


52  =  sin  E  (1) 

53  =  cos  E  sin  A 

The  sun  is  assumed  to  be  stationary  relative  to  the  earth  during  the 
entire  flight  of  the  projectile.  However,  the  field  of  view  of  the  yaw 
sonde  changes  as  the  projectile  follows  the  tangent  to  the  trajectory; 
hence,  the  sun  may  be  visible  at  one  stage  of  flight  and  not  at  another. 

The  orientation  of  the  projectile  is  given  by  unit  vector  V,  whose 
components  are 

V  =  cos  Q  cos  R 

V2  =  sin  Q  (2) 

=  cos  Q  sin  R 

where  R  is  the  range  azimuth  and  Q  is  the  quadrant  elevation.  This 
unit  vector  is  called  the  velocity,  but  it  really  is  the  vector  along  the 
projectile  axis,  which  is  here  assumed  to  be  tangent  to  the  trajectory, 
i.e.,  the  projectile  flies  without  yaw.  This  is  a  slight  limitation  that 
will  be  discussed  later.  The  sun's  visibility  was  calculated  for  three 
segments  of  the  trajectory,  by  setting  Q  equal  to  QE,  zero,  and  -QE, 
representing  launch,  peak,  and  impact,  respectively.  This  assumes 
a  parabolic  trajectory,  where  the  descending  portion  is  a  mirror  image 
of  the  ascending  portion.  This  assumption  is  valid  in  this  project,  be¬ 
cause  the  difference  between  actual  anti  vacuum  trajectories  is  within 
the  range  of  error  in  locating  the  visibil.ty  limits. 
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The  field  of  view  of  the  yaw  sonde  is  a  cone  with  a  90*  apex  angle, 
the  axis  of  the  cone  being  perpendicular  to  the  axis  of  the  projectile. 
The  angle  y-'  is  measured  from  the  nose  of  the  projectile  to  the  sun 
vector  drawn  from  the  yaw  sonde  position  on  the  projectile's  axis. 

For  the  yaw  sonde  to  see  the  sun  f  must  lie  between  45*  and  135*. 

From  the  above  descriptions  of  S,  V,  and  ^  ,  it  may  be  seen  that 
the  angle  is  given  by  the  scalar  product. 

cos  ^  =  S  V  =  S2  VA  +  S2  V2  +  S3  V3 

Since  ^  is  known  to  be  45*  and  135*  along  the  two  edges  of  the  field 
of  view,  it  will  be  practicable  to  use  it  as  an  independent  parameter 
for  finding  the  sun's  position  (A,  E). 

We  note  that  E  has  a  known  range,  hence  it  is  assigned  a  set  of 
values  in  2°steps,  from  0*to  about  80*.  This  leaves  A  as  the  remain¬ 
ing  unknown  angle.  From  Equation  1 

S2  =  sin  E 

and 

A  =  arctan  (S3/Sj)  (4) 

Thus  A  can  be  calculated  if  and  S3  can  be  found,  and  this  can  be 
done  by  using  the  directional  cosine  law  for  S, 

1  =  S^2  +  S22  +  S32  (5) 


in  connection  with  Equation  3.  The  details  of  the  algebra  involved  in 
solving  for  and  S3  sure  given  in  Appendix  A  and  in  the  listing  for 
the  computer  program  called  YSFIRE,  in  Appendix  B.  The  above 
calculations  were  made  for  each  of  the  three  firing  sites  in  Table  1, 
at  launch,  peak,  and  impact,  for  30*  ,  45*  ,  and  69*  QE. 
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TABLE  1 


Geographic  details  of  firing  sites 


Site 

Latitude 

(North) 

Longitude 

(West) 

DLG  mange 
(hour)  Azimuth 

China  Lake,  California 

35.75* 

117.67’ 

-.155 

.148.0 

Wallops  Island,  Virginia 

37.83* 

75.48* 

+.032 

128.0 

Yuma,  Arizona 

32.83* 

114.25* 

-.382 

90.0 

The  punched  cards  from  YSFIRE  containing  the  (A.  E)  points  are  sorted 
and  those  points  discarded  that  fall  outside  the  solar  grid  which  will  be 
described  in  the  next  section. 

These  results  tell  where  the  sun  has  to  be  for  yaw  sonde  sightings. 
The  next  task  is  to  find  when  the  sun  will  be  at  these  positions. 

Solar  Grid  Calculations 

This  second  line  of  computation,  to  find  when  the  sun  is  in  certain 
positions,  is  based  on  a  new  type  of  graph,  called  a  "solar  grid.  "  The 
technique  used  in  Reference  3  was  extended  to  combine  on  one  graph 
the  sun's  azimuth  and  elevation  as  functions  of  declination  and  hour 
angles  for  each  latitude  under  consideration,  as  Figures  2,  3,  and  4 
show.  An  array  of  (A,  E)  points  as  function!  of  D  and  H  was  generated 
and  plotted  by  the  computer  program  named  SOLAR,  using  the  equa¬ 
tions 

sin  E  =  sin  L  sin  D  +  cos  L  cos  D  cos  H  (6) 

sin  A  =  -cos  D  sill  Il/cos  E  (7) 

cos  A  =  (sin  D  -  sin  E  sin  L)/(cos  E  cos  L)  (8) 

Each  line  of  constant  declination  (except  tne  solstices)  represents  two 
dates,  because  the  sun  makes  a  round  trip  across  this  grid  each  year, 
as  Table  2  indicates. 
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The  compute? 
obtained  from  *h 
3, and  4  show,  ft  . 
and  elevation  on  , -» 
translate  the  sun h. 
(D,  H)  numbers. 


gram  SOLAR  also  plots  on  this  grid  the  limit  lines 
FIRE  program  in  the  previous  step,  as  Figures  2, 
ow  evident  why  it  was  expedient  to  combine  azimuth 
t  graph,  because  this  makes  it  possible  to  graphically 
limiting  angular  (A,  E)  positions  into  date  ai  d  hour 
This  is  done  by  hand,  reading  off  the  intersections  of 


each  limit  line  with  successive  date  and  hour  lines  over  a  complete  year 
cycle.  These  (D,  H)  number  s  are  punched  on  cards  for  the  final  plotting 
of  the  year  chart. 


Many  of  the  limit  lines  cross  the  grid  from  winter  to  summer  sols¬ 
tice.  For  such  lines,  a  standard  set  of  22  dates  (at  every  5*  of  declin¬ 
ation)  was  established,  to  expedite  the  transition  to  (D,  H)  numbers, 
starting  with  1  January  at  22.  0  S  declination.  This  type  of  limit  line 
transforms  to  lines  that  cross  the  year  chart  from  January  to  December, 
e.g. ,  the  impact  limit  for  45*  QE  in  Figure  5. 


Other  limit  lines  do  not  cross  the  grid,  but  enter  and  leave  at  either 
solstice.  These  lines  transform  to  closed  loops  on  the  year  chart,  e.g., 
the  launch  limit  fer  60*  QE  in  Figure  5.  Oi  the  same  figure,  the  limits 
for  peak  and  30*  impact  also  form  loops  which  are  interrupted  at  the 
year's  end. 


Year  Chart  Display 

The  actual  ploiting  of  the  (D,  H)  numbers  is  done  by  the  computer 
program  YSYEAR.  Special  calendar  year  graph  paper  was  obtained 
expressly  for  this  part  of  the  project.  *  The  dates  D  were  changed  to 
day  number  by  reference  to  Table  3,  to  fit  the  dates  into  the  366  divi¬ 
sions  on  the  long  axis  of  the  graph  paper. 


To  produce  Civil  Standard  Time,  the  local  hours  H  are  corrected 
for  Meridian  Passage  (MP  or  EMP)  and  longitude  in  the  YSYEAR  pro¬ 
gram  according  to  the  equation 


CVT  =  H  -  EMP  +  DLG  (9) 


*The  special  graph  paper  was  taped  to  the  paper  on  the  drum  of  the 
Calcomp  Plotter,  to  draw  ihe  graph  directly  on  the  year  paper. 
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*  Id  leap  jnm,  alter  February  28,  add  1  to  the  tabulated  number. 


i 
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The  MP  figures  are  obtained  from  Table  2  and  converted  to  hours  be¬ 
fore  reading  into  the  program,  thus 


EMP  =  (MP  -  12)/60 

The  MP  is  the  local  mean  time  at  which  the  sun  crosses  the  meridian 
at  any  given  site.  Throughout  the  year  it  varies  from  +15  to  -15  min¬ 
utes. 

The  correction  for  longitude,  DLG,  converts  local  time  to  civil 
time,  i.e. ,  the  clock  time  in  the  particular  time  zone  where  the  fir¬ 
ing  site  is  situated.  The  limit  lines  plotted  by  YSYEAit  are  shown  in 
Figures  5,  6,  and  7. 

Trajectory  Heights  at  Different  Terminal  Slopes 

A.  preliminary  review  of  the  windows  'or  the  different  quadrant 
elevations  in  Figures  5,  6,  and  7  showed  that  the  60’ quadrant  eleva¬ 
tions  are  the  most  restricted,  but  they  are  also  the  trajectories  that 
are  of  the  greatest  interest.  Since  the  descending  branch  of  a  trajec¬ 
tory  seems  to  be  less  important  than  the  ascending  it  may  be  worth 
estimating  how  much  of  a  trajectory  is  recorded  if  the  projectile  is 
fired  at  some  time  point  in  the  obscured  region  for  the  actual  quadrant 
elevation.  That  is,  how  much  of  a  60’  trajectory  will  be  recorded  if 
it  is  fired  at  a  point  on  the  impact  limit  for,  say,  a  45 ’trajectory. 
This  may  widen  somewhat  the  windows  for  the  steeper  trajectories, 
if  part  of  the  descending  portion  can  be  disregarded.  A  few  calcula¬ 
tions  have  been  made  to  investigate  this  possibility,  assuming  sym¬ 
metrical  parabolic  trajectories. 

From  elementary  mechanics,  che  equation  for  a  parabolic  tra¬ 
jectory  with  origin  at  the  launch  point,  is 

2  2 

y  =  x  tan  Q  -  gx  /(«?v  cm  ’  Q) 


The  rarge  is 


r  =  (2v  '  /g)  sin  Q  cos  Q 


and  the  maximum  height  is 

h  =  (v2/2g)  sin2Q 
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Through  suitable  manipulations,  the  velocity  is  eliminated  and  dimen 
sionless  height  and  range  coordinates  are  introduced,  giving 

7/ h  =  (4/r)  (x  -  x2/r) 


The  problem  now  is  to  find  the  height  at  which  a  given  trajectory 
has  a  smaller  slope  than  its  initial  quadrant  elevation.  This  is  done 
by  difierentiating  the  above  equation,  which  gives 

y*  /h  =  (4/r)  (1  -  2x/r) 

where  y'  is  the  slope  at  height  h.  Eliminating  the  dimensionless  range, 
we  find 


y/h  =  1  -  (ry'  /4h)2  =  1  -  (y*  cot  Q)2 

where  r/4h  =  cot  Q  and  y*  is  the  arbitrarily  chosen  smaller  slope. 

Table  4  shows  the  results  for  three  simple  cases  encountered  in 
this  project. 


TABLE  4 

Heights  at  early  cutoffs 


Qirdrant  Elevations 

y1 

y/h 

60’ 

45’ 

2/3 

60* 

30* 

8/9 

45* 

30* 

2/3 
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A  proje<  tile  fired  at  a  time  point  on  the  impact  limit  curve  for  y* 
will  go  beyond  the  peak  to  the  above  fraction  of  h  before  the  sun  disap¬ 
pears  from  the  yaw  sonde.  Conversely,  a  projectile  fired  at  a  point 
on  the  launrn  limit  curve  for  y*  will  not  see  the  sun  until  it  reacnes 
the  above  height  on  the  ascending  segment  .of  the  curve.  All  projectiles 
fired  at  time  points  on  the  peak  curves  will  either  cut  off  there,  if  in 
the  terminal  regio..,  as  in  Figures  5  and  6  and  the  upper  part  of  Fig¬ 
ure  7,  or  w  11  start  only  at  the  peak,  as  in  the  launch  region  of  Figure 
7. 


RESULTS 

The  results  obtained  in  this  study  are  the  solar  grids  and  year 
charts  for  each  of  the  three  sites,  Figures  2  through  7  inclusive. 
Table  4  shows  the  height  at  which  early  cutoffs  in  steep  trajectories 
occur. 


DISCUSSION 


Limit  Lines 

The  procedure  in  YSFIRE  for  calculating  the  limit  points  (A,  E) 
leaves  out  the  positive  identification  of  the  visible  side  of  tne  boundary. 
This  identification  is  fairly  evident  from  past  experience  and  from 
the  geometry  of  the  firing  range  situation,  hi  ambiguous  cases  the 
visible  side  was  identified  by  using  Rvalues  of  50'  and  130*  in 
Equation  3  and  noting  to  which  side  of  the  original  A  points  the  new 
A  points  fall  for  the  same  E  values. 

In  general,  the  limit  lines  may  have  minima  and  maxima  in  the 
solar  grid  region,  i.e. ,  there  may  at  some  places  be  two  different 
A  values  for  a  single  E  value  for  a  given  quadrant  elevation  line. 

This  possibility  is  anticipated  in  the  YSFIRE  program.  A  corollary 
to  this  is  the  requirement  that  all  limit  lines  enter  and  exit  the  grid 
somewhere;  they  dare  not  end  inside  the  grid. 

The  effect  of  the  actual  yawing  of  a  projectile  would  be  to  widen 
the  limit  lines  into  bands.  To  delineate  such  bands  for  a  given  yaw 
angle,  say  5°,  one  would  use  in  Equation  3  four  values  of^»  namely, 
45°  ±  5°  and  135°±  5°  .  This  would  produce  a  pair  of  A  values  for 
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each  E.  When  the  two  sides  of  a  band  are  transcribed  to  (D,  H) 
values  a  corresponding  band  would  be  produced  on  the  year  chart. 
But  the  width  of  any  such  band  is  not  easily  estimated.  In  practice, 
if  one  fires  too  close  to  a  limit  line,  periodic  breaks  will  occur  in 
the  yaw  sonde  record  when  the  angle  exceeds  its  limits  during 
projectile  precession  about  the  velocity  vector. 

Site  Parameters 


Two  site  parameters  that  affect  the  results  sire  range  azimuth  and 
latitude.  We  will  discuss  first  the  azimuth,  which  has  the  greatest 
effect.  This  parameter  has  nearly  its  maximum  possible  spread  of 
90*  in  this  study.  It  influences  the  A  values  in  the  (A,  E)  limit  points, 
as  can  be  seen  in  Figures  2,  3,  and  4,  where  the  latitude  spread  is 
only  5*  .  When  these  limit  lines  are  translated  to  tne  year  charts  the 
resulting  windows  present  very  different  appearances.  The  nearly 
north-south  azimuth  at  China  Lake  provides  the  widest  windows,  where 
as  the  east -west  azimuth  at  Yuma  provides  the  narrowest  windows. 

As  the  latitude  has  a  different  sort  of  effect  one  may  infer  that  north- 
south  azimuths  will  always  yield  wider  windows. 

If  new  firing  sites  are  contemplated  they  should  be  located  where 
nearly  north -south  range  azimuths  are  available,  since  the  rexge 
azimuth  at  any  particular  site  is  rather  permanently  and  narrowly 
limited  by  geographic  and  demographic  restrictions,  and  cannot  be 
changed  to  suit  the  convenience  of  ordnance  experimenters. 

It  is  also  interesting  to  note  that,  if  one  fires  a  projectile  in  the 
opposite  sense  along  a  given  range  line,  the  curves  on  the  solar  grids 
and  year  charts  will  be  the  same,  except  that  the  launch  and  impact 
labels  will  be  interchanged. 

The  other  site  parameter,  latitude,  has  less  influence  on  the  limit 
lines,  through  its  effect  on  the  solar  grid.  On  all  grids  the  equinox 
line  (O*  declination)  always  passes  through  the  sunrise  and  sunset 
lines  at  90*and  270*  azimuth,  respectively.  At  noon  it  passes  through 

an  elevation  equal  to  the  colatitude  (90*  -  L).  Thus,  all  grids  have 
the  same  width,  but  lower  maximum  heights  for  higher  latitudes. 

The  latitude  effect  for  a  given  range  azimuth  may  be  estimated  by 
imagining  that  the  limn  lines  for  one  latitude  are  drawn  on  the  grid 
for  another  latitude.  Because  the  slopes  and  shapes  of  the  limit  lines 
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themselves  differ,  each  quadrant  elevation  line  would  have  to  be 
studied  in  detail,  a  tedious  procedure  not  warranted  in  this  study. 

In  general,  it  may  be  inferred  that  change  in  latitude  will  change 
slightly  the  sizes  and  shapes,  but  net  the  topology,  of  the  yaw  3onde 
windows. 

Year  Charts 


The  limit  lines  on  the  year  chart  are  almost  symmetrical  if  re¬ 
flected  about  the  summer  or  winter  solstice,  and  they  would  be  ex¬ 
actly  symmetrical  if  it  were  not  for  the  MP  correction,  which  fluctu¬ 
ates  during  the  year.  The  MP  has  a  different  phase  from  the  sun  and 
has  two  different  periods  and  amplitudes.  The  maximum  variation 
is  from  +15  minutes  in  February  to  -15  in  November.  This  produces 
a  "vertical  shear"  in  the  lines,  between  3  November  and  8  February, 
on  which  days  the  sun  has  the  same  declination  (15  S). 

The  sunrise,  sunset,  and  peak  limits  are  common  to  all  trajec¬ 
tories,  regardless  of  their  initial  quadrant  elevation.  The  peak  limit 
is  the  same  as  that  of  a  flat  trajectory.  Starting  with  this  QE  of  0* 
and  proceeding  through  30*  ,  45*,  and  60*  quadrant  elevations,  one  can 
see  on  any  year  chart  successive  decreases  in  visibility.  The  steepest 
trajectory  has  the  smallest  window  at  any  site. 

The  windows  for  the  steeper  trajectories  can  be  extended  somewhat 
by  firing  at  time  points  outside  their  visibility  area,  if  one  is  content 
to  lose  part  of  either  the  ascending  or  descending  branch. 

How  to  Read  the  Year  Charts 


The  limit  lines  are  labelled  with  quadrant  elevation  symbols  identi¬ 
fied  in  the  legend  near  the  right  side  of  each  chart.  In  these  charts 
and  in  the  solar  grids,  the  letter  T  is  used  for  "Terminal"  in  place 
of  I  for  "Impact"  because  the  symbol  I  is  too  easily  mistaken  for  a 
number  1.  Lettering  in  certain  areas  states  which  segments  -.p  laiknch, 
peak,  or  impact  --  are  invisible  for  all  quadrant  elevations  surround¬ 
ing  the  lettering.  Other  segments  and  quadrant  elevations  are  visible 
in  such  areas.  On  the  Wallops  Island  chart,  the  launch  segments  are 
visible  above  the  lines  labelled  "Launch  Limits.  "  On  the  Yuma  chart, 
all  launches  are  visible  above  the  lower  group  of  lines,  and  the  impacts 
are  visible  below  the  upper  group,  e.  g.  at  Yuma  the  60°  launches 
(60  L)  are  visible  only  after  1235  hr,  and  the  60’  impacts,  60  T,  be¬ 
fore  1030  hr,  in  midsummer. 
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Individual  Sites 


China  Lake  offers  the  largest  windows.  All  launches  are  visible 
c.ll  year  at  any  time  cjf  day,  except  for  a  short  period  for  the  60* 
launch  (60  L)  Ground  noon  in  midsummer.  Impacts  areninvisible 
over  the  noon  .iour  all  year,  except  for  0*  and  30*  during  the  summer. 

Wallops  Island  presents  a  simpler  pattern.  All  launches  are  visi¬ 
ble  after  1330  in  the  summer  and  noun  in  the  winter.  Impacts  are 
visible  between  September  and  April  in  the  afternoons,  but  between 
April  and  September  the  peak  cuts  off  ai  1620  hours  and  the  steeper 
trajectories  cut  off  at  earlier  t;mes. 

Yuma  offers  the  smallest  windows;  here,  the  only  open  spot  for 
all  trajectories  is  from  0900  to  1400  hours  from  November  to  Feb¬ 
ruary.  The  peak  is  clear  between  0800  and  1500  hours  all  summer. 
The  launch  segments  for  the  other  quadrant  elevations  are  visible 
after  1030  or  1230  hours,  and  the  impacts  are  visible  before  1230 
or  1030  hours. 

These  charts  are  offered  as  a  rough  guide  for  selecting  suitable 
firing  times.  It  is  not  recommended  that  firings  be  made  at  times 
too  close  (say  15  minutes)  to  the  limit  lines,  for  reasons  mentioned 
earlier  in  the  Discussion. 


CONCLUSIONS 

It  is  practicable  to  plot  yfar  charts  of  yaw  sonde  windows  for  any 
given  site  and  range  azimuth.  The  windows  are  widest  for  north-south 
azimuths  (China  Lake)  and  narrowest  for  east -west  azimuths  (Yuma). 
The  steepest  trajectories  have  the  narrowest  windows.  At  the  three 
sites  studied,  it  is  possible  to  obtain  yaw  sonde  records  for  all  tra¬ 
jectories  at  some  time  during  the  year. 


EECOMMENDAT  IONS 

To  expedite  scheduling  of  firing  programs,  the  year  charts  sub¬ 
mitted  here  should  be  consulted. 
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When  new  yaw  sonde  firing  sites  are  chosen,  preference  should 
be  given  those  whose  range  azimuths  are  most  nearly  parallel  to  the 
local  meridian. 
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Fig  2  Solar  grid  for  Unina  Lake 


Wallops  Island 
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Fig  5  Year  chart  for  China  ^ake 


Year  chart  for  Yuma 


appendix  a 


Algebraic  details  of  the  solution  for  s;  and  s 
Elimination  of  S}  from  Equations  3  and  4  gives 


SlV  +  V32'  -  2S1V1<P  -  EV2,  +  (p  .  Ev/  v/(l  -  E2)  =  O 


The  solution  for  is 


V  (P  -  EV  ) 

s,  =  -i— - L. 

1  V**.* 


Then 


,  *  (t .  <vi2  -  w,)2  -  y3wr 1/2 

(V2  (P  -  EV  )  )c 


s3  =  (P  -  ev2  -  SjVj)  /V 


points  :ntheecromp:sasZi”olhciP.li“lS  °ne  °f  the  io“r 

me  compass,  special  cases  arise  because  either  S  or  S 

’  1  °3 


become  aero.  At  Yuma,  for  instance,  where  R  = 


90*  and  V  =  O. 

i 


s3  =  (p  -  ev2)  /v3 

and 
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Although  these  solutions  are  straightforward  for  hand  computa¬ 
tion,  many  branchings  must  be  anticipated  for  machine  computation, 
as  the  listing  for  YSFIRE  in  Appendix  B  Shows,  The  possibility  of 
two  A  values  for  a  single  E  value  is  also  provided  for,  by  including 
BAS1  and  BAS2  among  the  answers  when  punched. 

Because  of  the  way  the  program  is  arranged,  the  punched  cards 
are  a  mixture  of  different  stages  for  a  single  quadrant  elevation; 
however,  they  are  easily  sorted  out  manually  to  ppdduce  continuous 
limit  curves  on  the  solar  grid. 
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APPENDIX  B 


Program  Listings 


YSFIRE 

SOLAR 

YSYEAR 
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PROGRAM  (SF IKE (INPUT »T APF5= INPUT ♦OUTPUT. TAPE6=0UTPUT«PUNCH, 


program  to  calculate  opt  rwuw 


_ M  *  .  ^  r*  r  t  «v»«»  >  . 

tNtEOck  otn 


♦  *V (3*3) *5(3*3) *SITE (8) * ASW (3) *LSH ( 3) ,LSC ( J> *LSM < 3) ♦ LCM (3) , A (3) * 
- VBTT)  VA7T3TTP5TT)'  7E  ST44T «  AUT3*2T2T - - - 

iotr  rtjRw«r-tiHTr  - - —  - -  - 

101  FORMAT  (8A10i 

10?  FORMAT  ( 1?F6V0 ) - 

103  FORMAT  (1216) 

~tOR  FORMAT  ^/44X* 4 7H  OP T I MOMFTR IK  C  TTMFS ~FOR-Y AW  SONDE  TESTS  -  )~ 

105  FORMAT  (/21X«6HLAUNCH*5X«4HPEAK»4X*8HTERM1NAL/22X*2HAS*6X*2HES*6X 

-- ^»3MtMA»BX>^PCVTTI  -  - -  -  — '  -  - 

106  FORMAT  (  6X*A6*  5X »2F8 . 1 »F9.2 ♦ 15* 1 3) 


1UI  runnm  |  n  JAf  jn"al  f  3F  I  U  #  1  /  J 

108  FORMAT  <//35X * 59HCALCULAT ION  OF  SUN  AZIMUTH  FOR  GIVEN  PSI  AND  SUN 

'  +ELT  VATTON  V12H  FOR  O.F.~=TF5.1/T - - 

109  FORMAT  (  9X*A8«20X*30H  NO  SOLUTION,  IMAGINARY  ROOTS.  > 

111  FORMAT  <4X,3HPSI.16X,2HS2f8X*2HSP*8X*2HSM*8X,2HSl.HX.2HS3*  8X» 


112  FORMAT (/9X *  1 OHELEVAT ION  ,F5.1,J2X,  5HS?  =  ,F7.5,  44) 

-m  FORMA T  T4X.F4 . IT,  1  XVAftYSFl 0 .5 ,4F  10.1) - 

114  FORMAT  <12X,2F8.4,20X,2I4,4X,I4> 

-|-nr-pORMATtr30X ,T?MVEK IF  TCATION"  I 
♦CURVES  ON  YEAR  CHART  /) 


♦35X,57HPSI  MINIMUM  =  4b  DEG*  v  EARLIEST  TIMES  TO  FIRE 

117  FORMAT  (1H1,50X»18HINTERMEDIATE  DATES  *30X,6H0.E.  =,F5.I/ 

- *-35X*-<PWPST  MA XTMTIW  ^  T3S ~PF Or*  LATE S T  TIMES  TO~F IRF - 7 - 

7  6QX  *5HNST  =tI2) 


44X.F5.1) 


•ui^rotUf 


119  FORMAT  (20X * 6HR AS  1  = *F 7 . 1  *  1  OX , 6HBAS2  =,F7.1) 
C  30X  *  23H  AER08ALL 1ST  ICS  RRANCH  ) 


TTOXTArO. 


OATA  M0N/6HJAN  1*6HJAN  21*6HFER  8»6HFEB  23,6HMAR  8*6HMAR  ?1, 


uinnzmiJk'WJ  n*i  miii  me 


B  6HAUG  27* 6HSEP  10»6HSEP  23,6HUCT  6*6H0Cl  19*6HN0V  3*6HNOV  21* 
c  6  h  nrrc  ^  htje  c  tit  ~ 


:i»n. 


i. 


MC  =  1  TO  CALCULATE  AZIMUTH  FROM  GIVEN  PSI  AND  ELEVATION 


WRITE  46,1004)  JO 

WRITE  t6* I04T - 

READ  (5*101)  SITE 


iir>l  lilt  iiiTn:i:rj 


READ  (5,103)  MC,NCRV,KLT,NPCH 


WRITE  (6,101)  SITE 
IT  t  HCvtTQ  . 1  T~  GO  TO 
WRITE  (6,115) 

mtr-?? - 


rjlic.i  JuMCfi  l>f4Mlp14llVJniLf  I 


IF  (EOF (5) .NE.O  )  GO  TO  201 

- wa  *  of ( i )  - 


?ar  OO  9T  ="TT3 - 

OER  *  OE (L ) /RAON 


I  lilL'M  1  lilTM  I 


V (L  ,2)  =  SIN (OER  ) 
- 9  V(L ,3r  -  COSTQER  T 


WRITE  (6,108)  QE(1) 


INLE  =  14 


PS (2)  =  135. 


00  37  K*1,KLT 


ESR  =  ES (K ) /RAON 


WRITE  (6,112)  ES;K)  ,  S2»  ES(K) 


S3  =  0. 


ASP=  0.0 


JSP  =  PS(J) 


ft**  YSFTRF - TRttCE  — 


IT  TlTEtIiTTflN0"»H0VEQri'3tn  60 - TO  35 

IF  (I.E0.2.AN0.MQ.NE.60)  GO  TO  3S 


VA  =  V(L.l) 


VC  =  V (L.3) 

R  =  C0S(PS(J)/RA0N)-S2*VB 


0  =  1.0  -  S?*S2 

IF<RAZ.FQ. 90.0. OR. RAZ.EQ. 270.0)  GO  TO  211 
TFT R AT. ETT.”  OTO .t)R. PIT? .TO.  in 0  .  O T  GO  TO  712 


RT=  1.0  -  P 

I r  fvr.GT.O.O)  GO  TO  22 
WRITE  (6.109)  A7 ( I ) 

GO  TO“  3S  —  - 


k  i  =  qUH i  i  k  i 
S?  =  RV 

SP  =  RV*<1.0  ♦  RTF - 

SM  =  RV*  <1.0  -  RT) 

IF  <  RF.  EQ  .  0  .  OF  GO~T  (TTT 


S3M  =-<R  -  VA*SM)/VC 


IF  ( ASPR  . GE.0.0)  ASP  =  360.0  - 


27  SI  =“ST - - - 

S3  =  -  (R  -  VA*S1)/VC 


i  I  3J. 

IF  (AS. GE.0.0)  AS  =  360.  -  AS 


AQ(I.J.l)  =  AS 

- AQTI7J.2J  =  AT - - - 

GO  TO  29 

211  S3  =  -  R/VC 


J  I  —  W  “  3J" 

IF  (S1.GT.0.0)  GO  TO  210 

“WRITE  (6, Tim  A7(T) - 

GO  TO  IS 

~?TTJ  ST  =5QRTTSTT - 


c 


S3  =  0  -  S1*S1 


LW/411 


WRITE  16*109)  AZ(I) 


UU  I U  JD 

220  S3  =  SORT  (S3) _ _ 

AS2  =  ATAN2I-S3.S1)  «  RADN 
213  IF  IASI. GE. 0.0’  AS1  =  360.0  - 


WRITE  (6*1)3)  PS(J)*A2(I)«SZ*SP«SM*S1 »S3* ASPR* ASMR* AO  ( I  * J« 1 ) 

f  - - 

INLE  =  INLE  ♦  1 


=  AO ( I » J* 1 ) 


IF  (BAS1.LT. 60.0. OP. BAS1.GT. 300.0)  GO  TO  35 


WRITE  (6*119)  BAS1.8AS2 
IF  (NPCH.EQ.0)  GO  TO  35 


WRITE  (7*114)BAS2.ES(K)» I *MG, JSP 


35  CONTINUE 
^3~C0NTINUE  ~ 

INLE  =  INLE  ♦  2 


WRITE  (6.100) 


INLE  =  5 


WRITE  (6.100) 


t - VER1 FTCATTON  "OF  -SUN- 

C 


51  READ  (5.118)  DYN.OUN.0AS.BES.BSHA. BMP. LP.KTP.MQP.NSP 

- JF— =  TO - 

IF  (EOF (5) .NE.0)  JF  =  1 
- TT  fOF.EQ.n - GO  T0~T2 - 


8  WRITE  (6*105) 


tr - - - 

49  S<1.1)  S  COS (BES/RAUN) *COS  (BAS/RAOn) 


~ STT72T”3’5INfBES/RAON> - 

5(1.3)  =  COS (BES/RAON) *SIN(8AS/RADN) 

~c 

DO  S3  L=1 *1 

l)l>  *>2  1  =  1  *  J 

A ( I )  =  S(1*I) 

c 

~5?TTTT )  -  VTUTT  ) 

CALL  ANU TA «B *PSI (L) ) 

PSI(L)  =  RA0N*PSI(L) 

c 

53  CON.INuE 

t - CTVTL  TIffE — CALCULATED "FROM  LOCAL  HOUR  ANGLE  *  HER  11)  I  AN  PASSAGE* 

C  AND  DLG  FOR  EACH  DATE 


IF  (NST.EQ* 1 )  GO  TO  54 

- BCVT  -  BSHA  -EMPP4)  ♦  DLG - 

GO  TO  55 

59  BCVT  g~B5HA  -  EMP  ♦  DLG - 

55  LSHB  =  BCVT 

- LSCB  *  LSHB  *  T(n) - 

LSMB  =  BCVT  *  100  -  LSCB 
- LCMB  a  LSMB«60/I00 - 

- IF  (NS T7ETJ7DT  'D/N'  '-  "HON  CM) - 

WRITE  (6*106)  DVN*BAS*BES*BSHA»LSHB.LCMB 

- WRTTE  (6.10? T - TFSTT — LJ  .1=1.11 - 

INLE  =  INLE  ♦  3 

- rr-T rNUF.TT.641 - DITTO  T6 - 

WRITE  (6*100) 

- ,NLt  s-TT - 

1 6  TT  =  "*R  1  I 
GO  TO  51 

[  T7TC“=HNTT*T - 

i  IF  ,NC.GE.NCRV>  GO  TO  201 

[  - IF  (KIG.EU.2.AND.JF.EU.1) - GO  (DFOTJ - 

GO  TO  50 


201  STOP 
"  END 


SUBROUTINE  ANU( VI »V2r THETA) 

t - 

CANU  ANGLE  BETWEEN  TWO  UNIT  VECTORS 

•c - 

DIMENSION  VI <3  )»V2(3  ) 

■c - 

0P*V1U)  *V2<!)  *V1<2)  «V2<2>  ♦VI<3)  *V2I3) 


ULtKlIlj 


i  rL  v  r 

«5  RETURN 

- - - END - 


o  n  n  r>  n  r> 


PROGRAM  SOLAR ( INPUT .  T  A PES= INPUT. OUTPUT » 1  APF6=OUTPUT ) 

PLOTS  YAW  SONDE  LIMIT  LINES  ON  SOLAR  GRID  FOR  GIVEN  LATITUDE  ANO 

range  azimuth. 

GALLOPS  T SLANG  CURVES. 

U I MENTION  Hill ) .SNHIHl ) , CSH < P 1 > , A (} 2 » H? > . t ( 12.82) , DEC < 1 1 ) .HR ( 1 7) * 

♦  A  a  (PI)  «EF  m )  .SITE  (8)  .FS  (HO), AS  (HO).  AZ  ( 3>  ♦  JQ  (Hi  >  * 

2  KR(ll).KE<m.PO(?) 

C 

100  FORMAT  ( l?Fh.fi) 

101  FORMAT (/SSX.21H  LOCAL  HOUR  ANGLES  / 12X *  1  )  ( 7X . A2 ) I 

10?  FORMAT  ( 7X » 1  Ah  DECLINATION  =.  AS»4X.AHKPT  =  »I3.p/*4HKET  =  . 13.4X.4hNK 

♦  = ,  U.6X  »4HJ  =,I3,6X.3HJ0=.I3  /3X . 9HELLV AT  ION .  9F9.4/ 

♦  3X.9H  AZIMUTH  .  9F9.4Z) 

103  FORMA1  (  2X.HA10) 

104  FORMAT  (SX.7H  SlNt  H.1JF9.4/) 

105  FORMA!  (/TlX.S^h  SUN  ELEVATION  AND  AZIMUTH  VS  DECLINATION  AND  HOUR 

♦  ANGLE  /) 

107  FokmAT  <//?lX,n,??H  LIMIT  CURVES  PLOTTED  > 
m  FORMAT  (10X.PF1 0.4,31 3 .4X.A6) 

1 1 4  FORMAT  (2X.4F 10.4, 3X. 313) 

1 1 5  FORMAT  ( 4 1 S ) 

HA  FORMAT  (1?X.?F*.4«?0X»  n.?X,A2> 

119  F0RMA1  (/53X .?4hGR ID  POINTS  Al  RIGHT  F  NO  /PX.SpKET  =  » 7x . I  2 » 1 0 ( PX » 
9  I? 1 /8X , 7HFLFV •  , 1  IF 10.4/HX.7FAZIM.  .11F10.4) 

1004  FORMAT  (  1 H 1  //  45*  ,  IPHFNGINEFRING  SCIENCES  LABORATORY  /10X.A10. 
r  30  X  ,  P  1H  At  RfjRALL  1ST  ICS  BRANCH  //) 

C 

DATA  DEC/SH?3.4m,3HP0N  .3H15N  .3H1DN  .3H  SN  ,3H0.0  .3h  SS  , 

♦  3R 1  OS  .3H1SS  .3HP0S  .SH23.4S/ 

C 

DATA  HR/PH4  .?HM  ,?H6  ,?H7  ,?F«  »2H9  , ?H 1 0 ♦ ?H11 ♦ 2H i 2 » 2H1 3 , ?H 1 4 » 

♦  ?H|S.2H16-.?H17.?H1H.?HI9,2H20/ 

C. 

DATA  A/9H4«0.0/ 

DATA  c / 9 8 4 » 0 . 0 / 

DATA  AZ/4H  L , 4 mP  ♦ 4H  TZ 

DATA  hO/ftHj  C0NS.6HK  CONS/ 

C 

DATA  IV?3.4,?0.0, 15.0, 10. 0,5. u» 0.0, -5.0, -10. 0,-1 5.0. -20.0,-?  1.4/ 

c 

wADN  =  57.29S7P 
CALI  DATF  IJP* 

CALL.  PLOT  (.1,0, 5.0,-  I) 

1  WRITE  (6,1004)  JD 
WRTTE  (6.105) 

C 

READ  (S.103)  SITE 
READ  (5,100)  ATO.DLG.RAZ 
IF  (FOF (51  ,NE  .0.0)  00  TO  311 
READ  (5,115)  KW 
A  T (JR  =  ATU/RADN 
SAT  =  SIN ( AT UR ) 


Reproduced  from 
bes)  available  copy. 
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•gram 


.SOLAR 


TRACE 


.CDC  6400.  FTtt  V3«fl-£239  QEI=1._4 


COI  =  COS  (ATUR) _ _  _ _ 

C 

WRITE  (4.1u3)  SITE 

WRITE  I*. 101)  (HR  <L )  «L= J  » 1 7.P ) 

C  _  ...  .  .  _ 

CALL  SYMBOL  <0.7,6.5,S3H  SUN  ELEVATION  AND  AZIMUTH  VS  DECLINATION 

-  _ t_ AND  HOUR  ANGLE _ i5B»0..21jjQ»01_  ..... _  . _ _ 

C 

...CALL  SYMBOL  < 0.5,6,  1  »S I TE, HO* 0 . 17, 0. 0 )  .  . 

call  SYMBOL  <1.7,5.6.74HSH0WING  LIMIT  LINES  FOR  LAUNCH,  PEAK,  AND 
♦IMPACT  AT  0«E.  =  30,  45.  60  DEG.  ,74,0.14,0.0) 

C 

00  ?  K  =  l.« B 1  _  .  _ _ _  _ 

H  =  <57.0  ♦  K*3.0>  -  180.0 
SNH(K)  =  S IN ( H/P AON ) 

2  CSH  <K  )  =  COS  (H/RALIN) 

WRITF  (4.104)  (SNH(K) .K=l,8l,10) 

c 

CA.LL  A  A I S _ <JL*  0,0.0  ,.l  3HA/ 1 MUTH  <D£GJj13,  IP.  Gil).  Q  ±60^1^20^0  L  _ 

CALL  SYMBOL  < 0. 7 , 0 .4 , 1 7H0ECL IN AT  ION  (DFG ) « 1 7,0. 12,  0.0) 

CALL  PLOT  (0.0, 1.0, -3) 

EJS  =  -0.3? 

EJF  =  -0.1 
JE  =  1 

KBt=  1  _ 

KET=  1 
C 

00  517  M=1.81 

517  JQ(M)  =11  .. 

C 

.  .  _D0  .59  J  S..LUL  ■  -  —  -  _  _ _ 

C 

OR  =  r>13)/RAON  ..  __ 

SID  =  S  IN (OR ) 

CQO  =  C.OSIDR)  .  ... 

c 

C  SYMMETRY. OF  GRID  USED  TO  CALCUL  A  IE  .RJGHT  SIDE  FR.QM  LEFJ.S.UIL. 

C 

DO  4  K=l,81  _ 

IF  (K.GT.41)  GO  TO  57 

.  SNE  =_SAT*SID  ♦  CUT*CODaCSB(K)  . 

C5E  =  SORT (1,0  -  SNE*SNF> 

..  _  _  E < vL»Kl_  3  R AJ1N*  ASJP(SNF)  .  _  .  ..  _  .  _ _ _ 

E ( J, 42-K )  =  F ( J,K) 

C  .  .  ...... 

SNA  =-  C0D*SNH(K )/CSE 

CSA  =  (SID  -  'iNF*SAT>/<CSE*COIl  .  ...  _  .  ......  . 

A  ( J » K  )  =  RADN«ATA,'(?(SNA,CSA) 

_ IF  -iA10iiLL.iJ.0.QJ  _  A(0,K)  =  AijjK J_t  M.fl  .  _  _ 

A ( J.B2-K)  =  300.0  -  A ( J , M 
C 

57  CONTINUE 

4  CONTINUE  _ 

59  CONTINUE 
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o  o 


GRAM 


SOLAR 


COC  64UO  FTN  \)3.0-P239  OPT'l 


CALL  S04LF  <  A  *4H4  •  1 .12.0* 6Ql.0»2O«O  ♦  1  > . 

CALL  SCALE  <F*9m<**1.  4.0*  0.O*20.0»l> 

AN*3  =  0.0 
ANF  =  J2.0 

DO  61  J=l*ll 
00  60  K  =  t  *81 
A A (K  )  =  A  < J»K  ) 

EE  tK )  =  F  ( J.K  I 

IF  (K.EO.l)  60  TO  60 
IF  (K.GT.41)  GO  TO  55 
C 

C  J  CINSTANT 

C 

'll  IF  <F  ( J*i\  +  )  )  .LT.O.OO.Ow.FIJ.K  >  .GT.O.OUl  (  0  TO 
KMT  =  K  ♦  1 

SLJ  =  (FTJ.KPT)  -  £U«6Dl-in/(A(J»KHT>-A<j,KBT-l>> 
AjRT  =  A(j.KHI-l)  -  F < J.KHT-1 )/SLJ 
KET=  MR  -  KHT 
AJF  r  =  12.0  -  AJHf 
C. 

K  CONSTANT 


S4 


61 


IF  (F(J.M  .LT.O.OO.OK.E(J*1«K).GT.O.(K')  GO  TO  60 
JF  =  J 
JO  IK)  =  JF 
JO<H?-,<)  =  JF 

IF  (K.GT.KHT  .Aivn.  J.tU.  *  1  )  GO  TO  60 

SLK  =  <  F  ( JF  +  1 « K  >  -  t  (JE  ,K  )  )/ (A  <  JE*1  *K  >  -  A(JF,K>) 
ANP  =  A(JF.F)  -  t ( JF  »K 1 /SLK 
ANF  =  12.0  -  ALh 


c 


KJ  T  =  6 

t\FT  =  -  KMf 

IF  (J.Fg.  ID  GO  TO  60 
CALL  BLOT  <A(J.i\r)T)tE(J*KbT>»3> 
CALL  PLOT  (  AN  H  •o.0*R) 

IF  (J, FQ.il)  GO  TO  60 
CALL  FLO  T  (  A  <  J  *  K  E  T)»F  T  J  *  K  E  T  )  *'5) 
CALL  PLOT  ( ANE *0.0*2) 
ss  CONTIMIJF 
60  COM!  INI  IF 


Reproduced  from 
besl  available  ci-nv 


NK  =  KE 1  -  ^hT  +  I 


IF (FE(KHT).lT. 0.001)  GO  TO  6? 
CALL  PLOT  <aa<KPT)  . FF <  6  H  T  > ,3) 

CALL  PL-  T  (AJhl  *0.0.2) 

C 

C  WkITF  DtCLlNA I  IONS 

C 

62  CONTINUE 

c 

6  IF  (J.LQ.l)  GO  TO  T 

IF  <MuO ( J.2) .NF .0)  GO  10  62 
A JS  =  A JbT  -  0.1  I 
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o  o  n 


SOLA* 


Tk  ACE 


CDC  64UO  FTN  V3.0-P239  OPI=l 


CAUL  bYMHOI  (AJS*£JS»OEC(J).  1*Q.I2*  53. OJ 
IF  < J.F0.4.0R.J.EQ.6.0K.J.Ef ,P>  GO  TO  j 
GO  TO  54 

.1  XJ  -  A  A  ( 4 1 )  -  0.14 
Y J  =  EEtM)  -  .Oh 

CALL  SY4M0L  (X J.YJ.OFC'  1 » »4, 0. 12» 0.0 > 

58  CONTINUE 

IF  IJ.NF.l 1)  GO  TO  h4 
5JX  =  AU4T  -  O.hk 
SJY  =  -  ' .49 

CALL  bYi'MOl  <S.JX«SJY.OEC(J>  »5»0,I2*  53, 0> 
C. 

64  CALL  L INF  (AA(KhT) .£E<KrtT) ,NK.i,0»U»0) 

C 

IF (EF (KFT) .LT.O.OOi  >  60  TO  63 


CALL  SYMBOL  <AJF,tJF»OEC(J)»b»0.12*“b3.0) 
C 

5  CONTINUE 

IF  (J.NF.ll)  GO  10  61 
S J  x  =  A JFT 

CAI  L  SYMBOL  ISJX.EJF»0EC(J>»5»U.1?«-C»3.0) 
C 

61  CONTINUE 
C 

L  =  1 
C 

K6T  =  KHI1) 

NET  =  KU1) 

00  9  K-KHT.KE1 

CALL  LINE  (A(l.K)  itl  i  .  K )  *  10  IK  ).!,<», 9.0) 
IF  IK.EO.l)  GO  TO  10 
M  =  K  -  1 

IF  (MUO (  M.bl.NF.O)  Gu  TO  9 
10  CONTINUE' 


Mb’  n  L  HOUR  AMGLtS 

IF  (  K.NF.41)  GO  TO  is 
XU  =  A  < i * 4  l ) 

YU  =  E  1 1  .**  I  »  ♦  ,1 
XL  =  A(ll,4l) 

YL  =  E  (11 .4  1  )  -  .1 
ABU  =  XU  -  .Oh 
EPU  =  YU  ♦  .01 
ABL  =  XL  -  .Of 
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gram 


SOLA* 


TkACF 


COC  6400  FTN  v3.0-P239  <JPT  =  1  0 


C 


c 


c 


c 


c 

c 


c, 

c 

c 

c 


FHL  =  YL  -  .11 
GO  TO  I  / 

IS  SL^'I  =  <F(?.K)  -  t(l.K))/<M?.M  -  A(1.K>) 

SLPL  =  OKlliM  -  t(10*K))/(A(ll*K)  -  »(lo »K)) 
IF  iSLPL.FO.O.O)  GO  TO  9 
XL  =  A ( i 1 . K )  -  ,1/SLPL 
yl  =  em.K)  -  ,i 


IF  (  K.GT.41 )  GO  TO  )  1 
IF  (SLPU.FO.O.O)  GO  lu  9 
XL!  =  A(1,K>  -  .1 
TU  =  E(l.K)  -  . 1  *SLPU 
API)  =  XI)  -  .1)4 
FRU  =YU  ♦  .02 

IF  (K.iME.31)  GO  TO  14 
aha  =  AHI)  -  0.4 
EGA  =  tHIJ  -  n.f, 

CALL  SYMBOL  (  a^A,  f  ha  .  1 7HL0C  Al_  HOUR  AMC-LFS  •  1  V ,  0 . 1  2 .  4b .  0 ) 
lu  AfiL  =  *L 

FPL  =  Yl.  -  .1? 

C-0  TO  17 

1 1  IF  (SLPU.EU.O.O)  60  TO  9 
Xu  =  A  (  1  ,M  ♦  .  j 
YU  =  F  < 1 «K )  ♦  . 1*9LPH 
ARu  s  XU  -  n.n7 
EPII  =  Y'J  +  0.02 
APL  =  XL  -  .IU 
EHL  =  yl  -  .13 

17  COMT I  NOE 

L  =  I  ♦  K/s 
CALL  SYMROt.  (  ARU.FhU.HR(L). 2.0. 12.0.0) 

CALL  PLOT  (XU.YU.3) 

CALL  PLOT  (Ail ,K> ,F  < I .«> .2) 

IF  (K.LT.KP(U)  )  GO  TO  9 
IF  (K.GT.KE  (ID)  GO  TO  4 
CALL  PLOT  (A(ll»K)»E(ll»K)»3) 

CALL  PLOT  (XL.YL.2) 

CALL  bTMBOL  (  ARL » FUL . HK (L ) . 2 . 0 .  12 . V  .  0 > 

4  CONTlNUr 
42  CONTINUE 

CALL  AXISi  D.OtO.O*  ISHFLF  VAT  ION  (DEG)t-]S«4.P.'?0.0.  0.0.20. 0> 

CALL  AXIS  <0.0, O.L  .  IShFLfcVATION  (DEG),  IS**..  0.90.0.  0.0.20.0) 

IF  (Kk.FO.O)  GO  TO  311 

LIMIT  CURVES  PLOTTED.  IF  M<.GT.2FkO 


IK  =  0 

14  M  =  j 

20  Kf  AO  ( S  *  1 1  rt )  AS(M) .FS<iD  . IIEMP.MUP 
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IGrtAM 


SOL  Art 


TbACF 


LHC  #i4U0  FTN  V3.0-P?3V  OPT=l 


r 


C 

C 

C 


r 

C 


c 


IF  (£0FI5) .NF.O.OTGO  TO  <!1 
I  =  I  TEMP 
rIGF  =  MOP 
hi  =  M  *| 

SO  TO  ?G 

21  J  =  M  -  1 
CALL  SCALF 

CALL  SCALt  LES.J.l.  u.O ♦  U. 0.20. 0.1) 


CALL  5YM«OL  <»'A.wY.MOE,2«U.l?.u.O> 

22  CALL  SYMBOL  (W«,nY,AZ(n  .A  .  U  .  1?.  ('  .0  ) 

CO  TO  2H 

<?4  CALL  SYMBOL  (WX»ivY«AZ<I).  4,0.i?.0.'J> 

?8  IP  =  IK  ♦  l 

fF  (IK.Lr.sw)  GO  10  19 

CALL  SYMBOL  JOilOHL  =  LAUNCH* 10. 0. 12.0*0) 

CALL  SYMBOL  <G.SK*,:|.Ub*  BHP  =  BE  A*  ,  n,n.  IP.  u.O) 
CALL  symbol  <S,4  .-.2o,12hT  =  TEPMINAL.12.0.12.G.0) 

v/BTTE  (ft.in.7)  f  k 

in  stop 
L\P 
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I 


PROGRAM  YSYEAR( INPUT. TAPE5= INPUT .OUTPUT »TAPE&=ODTPUT> 

r 

C  PkOGRAM  PLOTS  YA*  SONDE  alNDOWb  POP  WHOLE  YEA« 

C 

DIMENSION  Dyn(50> ,CVT(SO) ,hR(14) ,QE (3) ,EMP(?2) , DUN (22) 

♦,M0n(22> .STTF(R) ,NSY(6» 

C 

C  CHINA  lake  CURVES. 

C 

101  FORMAT  <4  IS) 

102  FOPMflf  U2F^.O) 

103  FOPMAT  (PAHO 

104  FORMAT  (  //  1 7  X ♦ 13 ♦ 15H  CURVES  PLOTTED) 

118  FORMAT  (  An  » r6 •  0 ,4FH .4 ,4 1  ^  ) 

IIP  FORMAT  (//43X.36H  year  CHART  OF  YAw  SOnCE  wlNUOwS  AT  /) 

1  Pm  FORMAT  ( 2 < 1 0* ,4HQATF *  1 1 UX * A6 ) / I  OX  * 4H  CV I  *  1 1 (F 1 0 •? > / > / > 

100a  FORMAT  (  1 H 1  //  45X  .  32HFNGINEFKING  SCIFnCES  LABORATORY  /10X.A10. 
C  3(>a.  2  1h  AEPOPALLIST ICS  BRANCH  ) 

C 

DATA  PR/2H  6. PH  /,2b  8 .  PH  9.2Ht',2Hll,3P12.2HH,2Hla,2'-<>B,2H16. 
♦2P17,2hl8,?H19/ 

C 

DATA  EMP/o. 0No7, 0. Irf33.. ).? 333. 0. 2167. 0.1S-13.0.  1)67,0.066  7,0.0. 

A-!).  oS.-o. 0667,0.0  333.0. 10,0.02  i3.0. 0333,-0 .05.-0 •  1  3  <3,-o •£(> , 

P-0. 2S.-0. 26-7,-0.2333, -0.0333. 0.0333/ 

C 

DATA  dun/  1.0, 

,  2  1  •  0  »  ^9 . 0 ,54 . 4 ,68. 0 , 8  t  .  0  »  V3  »  0  »  1 0  7  •  0 , 1  22  •  0  *  14  1  •  0 , 1  73.  (I  , 

C  20S.0* 225. ().P4rt.0»2Pa.u,?67.0, PRO. 0,293.0, 3 08.0*3?6.0»3S6.0. 
f)  366.0/ 

C 

DATA  M0N/6HjAN  1»6HjAn  21.6HFEH  P.6HFEH  23*6HMAR  P»hhi*lAP  21, 

A  tHAPR  2 , 6h APR  16.6HMAY  1 , 6PM A Y  20,bHJUN  21.6PJLL  ?3*6HAIJG  12. 

P  6HAIJG  27.6HSEP  10.6HSEP  23,6P0CT  6.6HQCT  19.6PN0v  3.NHN0V  21, 

C  6HDEC  2 1 »6HDFC  31/ 

C 

DATA  NSY/03,04,06. 00.62.S1/ 

C 

CALL  DATE(JD) 

CALL  PLOT  (0.0, 0.0, -3) 

CALL  SYMBOL  (0. 0,7. 2, JD, 10.0. 14,0.0) 

CALL  SYMBOL  (2. 5,7.2, 3SHYAW  SONDE  WINDOWS  FOR  WHOLE  YF.AR  ,35,0.  IP 
♦ ,  0  •  0 ) 

7  READ  (5,103)  SITE 

IF  (EOF (5) .NF .0.0)  GO  fO  49 

READ  (S.102)  ATU.DLG.RA7 

CALL  SYMBOL  (0.3.S.7,  S I TE ,80 . 0 . 14, 0.0 ) 


C 


c 

LU 

1 

MU  30 

c 

LQ 

2 

MQ  45 

c 

LR 

3 

MO  60 

c 

LQ 

4 

peak 

c 

LU 

5 

SUNSET 

c 

LQ 

6 

SUNPISE 

c 
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o  r>  r>  on  o  nononono 


KT G  =  1  FOR  POSITIVE  SLOPE  (LAUNCH) 

KTG  -  ?  FOR  NEGATIVE  SLOPE  (IMPACT) 

KTG  =  3  FOR  SUNSET  LINE 

NST  =  0  FOP  STANHARO  DATES  (22  OF  THEM) 

NST  =  1  FOR  INTERMEDIATE  OATES 
NCRV  =  NUMBER  OF  CURVES  TO  BE  DRAWN 

READ  (5*101)  NCRV 
WRITE  (6*)004)  JO 
WRITE  (6*119) 

WRITE  (6* 1 03)  SITE 
NC  =  0 
9  1  =  1 

11  READ  (5*118)  AFD*DYN(1)*BAS»BE3*SHA»BMP*LIP*KTP*M0P,NSP 
IF  (EOF(5)  .NE.0.0)  GO  TO  IP. 

LQ  =  LIP 
KTG  =  KTP 
MQ  =  MQP 
NST  =  NSP 

IF  (NST.EO.l)  GO  TO  21 
CVT ( I )  =  SHA  -  FMP(I)  ♦  DLG 
DYN(I)  =  DUN ( I ) 

GO  TO  22 

21  CVT(I)  =  SHA  -  BMP  ♦  OLG 

22  I  =  I  ♦  1 
GO  TO  1) 

12  N  =  I  -  1 

CALL  SCALF  (DYN*N*l*l«j. 0*0. 0*36. 6.1) 

CALL  SCALE  (CVT.N.1,7,0*  6.0,1.818,1) 

CALL  LINE  (DYN.CVT *N* 1 *NSY (LQ> *0. 10*0) 

IF  (LQ.FQ.6)  GO  fO  4<j 
IF  (KTG.NE.3)  GO  TO  41 

CALL  SYMBOL  ( 8 . 70 * 6 . IS , 6HSUNSF T * 6 , 0 *  1 0 ♦- 16. 0 ) 

GO  TO  41 

40  CALL  SYMBOL  ( 8 . 30 , 0 .5  7  *  7HSUNR  I SE  *  7 , 0 . 1 0 , 1  7  .  l) )  . 

41  NC  =  NC  ♦  1 
IF  (NC.LT.nCkV)  GO  TO  9 

DRAW  HOIJP  L  INFS 
DO  4b  K=1 » 1 3 
Yh  =  U.55  *  (K  -  1) 

YL  =  »H  ♦  .0? 

YR  =  YH  -  .07 

CALL  SYMBOL  ( 0 . 05 * YL .HR (K ) ,2, 0 . 12. 0 . 0 ) 

CALL  PLOT  (0.0.YH.3) 

CALL  PLOT  (10.0. YH. 2) 

CALL  SYMBOL <1 0.04, YH,HR(K) ,2,0*12.0.0) 

45  CONTINUE 

C 

C  DRAW  VFkIICAL  AXES 


40 


c 


c 


c 


C^LL  PLOT  ( 10.0*6.65*3) 

CALI  PLOT  <10.0.0.0,2) 

CALL  SYMHOL  (-0,05,2.5, 19HC1VIL  STANDARD  TIMF » 19*0. 12*90.0) 
CALL  PLOT  <  0.0.6.85.3) 

CALL  PLOT  <f. 0,0.0.21 

CALL  SYMBOL  <3.«,5.?,20HaLL  SEGMENTS  V ISIHLF ,20, 0. 1a ,0 .0) 
CALL  SYMBOL  (  j.B*0.6.20HALL  SEGMENTS  \l  TSIbLF  *  2G  *0.  14  ,0 .0 ) 
CALL  SYMHOL  (  7 . 1  ,  2 . 9 , 22h30T  AND  PEAK  INV  IS  l  BLF  ,  2?  , 0  .  1  A,  0 .0 ) 
CALL  SYMROL  (  <) .  7 »  3  •  0  »  22HPE AK  ANO  30T  IN  V I S I  RLE  *22 . 0 . 14 , 0 .0) 
CALL  SYMBOL  <3.95,3.5,  13H60L  INV  I  SI  »Lt  *  1  3  *  0 . 14 . 0  .0 ) 

CALL  SYMBOL  (3. 5,1. 8,  1  7HIMPACTS  INv  ISlRLF.  , 1  7 , 0 . 14,  U  .  0 ) 


C«LL  SYMHOL  (9 
CALL  SYMBOL  (9 
CALL  SYMHOL  (9 

call  symbol  <y 

CALL  SYMHOL  (9 
“oBTTF  ( 6  «  1  0  A ) 

A 9  STOP 
END 


20, 

3. 

.90  , 

7HSYM 

20* 

3. 

.715. 

7H 

n 

20* 

3. 

.577, 

7H 

3 

2U. 

3, 

.439* 

7h 

4 

r' 'J  , 

3, 

.301, 

7H 

b 

NC 


QL. 7,0.10*0.0) 
0,7,0.10*0.0) 
30*7,0. 10»U.0) 
45,7,0.10,0.0) 
60*7*0.10*0.0) 


fieProduced~froin 

^£L»vailable  r»n. 


